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ainhal Posi tinnmp Systems 

The present invention relates to developments intended to increase lire accuracy 
obtainable trom global navigation satellite systems (GNSS). 

At present mere are two publicly available GPS systems, known as NAVSTAR, 
5 owned by the USA, and GUONASS owned by the Russian Federation. These have 
been in existence for around two decades, but in the near future it is hoped that the 
European regional augmentation of GPS will start to provide its services, followed 
within a few years by a European system under the name of GALILEO. 

The existing systems have been progressively refined so that using a differential phase 
1 o implementation a location* accuracy of less man 2 cm can potentially be achieved 
ever a baseline of 1000 tan, hut with a cost in computation and in the time taken to 
determine the location. Real time or near teal time measurements have a 
correspondingly lower resolution, and at present the requirements for high precision 
mean mat additional augmentations are necessarily employed to supplement the 
15 GNSS information. Furthermore, these could include a receiver taking measurements 
from many satellites, up to all those visible to it whereby to calculate an over- 
determined position solution and rejecting inconsistent data to improve the accuracy 
of the position solution. Such a system may use data from more man one constellation 
of GNSS satellites, GPS and GLONASS. However, a significant error arises from the 
20 inability to tnodel accurately the delay to the GNSS signals caused by the troposphere. 

One of the largest identified sources of error lies in the effect that the troposphere has 
on signals propagating therethrough. The troposphere introduces ray bending and a 
signal delay which is influenced by a number of meteorological factors, but 
particularly water content 

25 Traditionally, fhis delay h,s been 

models. TTris weather data is too large to b* seat over commnmcatioii systems that 
are available to mobile users. 
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However, recently an attempt has? been made to introduce numerical weather 
prediction (NWP) itxto the GPS process, whereby local estimates of tropospheric 
delay T^^fie" pro^d«L ™ T^'fiSntiom was to use high resolution NWP for 
deteimining the tropospheric zenith delays for the estimated position of a GPS 
5 receiver, instead of using a global delay model, which delays can subsequently be 
used in the positioning process. 

However, even if the zenith delay is known, there will be a need for a user to be able 
to correct for other directions if increased accuracy of measurement is desired. 

The present invention provides a method of enhancing the accuracy of the GNSS 
.1 * . • • 
1.0 — posifioiLsolution.by.coTnm dam that enhaoce$ its -capability to~ 

model the tropospheric error. The method comprises the steps of generating a first set 
of less accurate tropospheric delay values from a first model which is known per se 3 
generating a second set of more accurate tropospheric delay values from a second 
model based on meteorological information, and developing a set of modified said 
15 parameters for use with said first model so that it thai provides a set Of delay values 
substantially approximating the second set, and expressing the modification to said 
first parameters as a set of correction parameters. In an embodiment of the invention, 
these parameters are differences with those obtained from the first model and are 
transmitted to the receiver, 

20 In one embodiment the first model is based on one or mote non-meteorologicai 
parameters). These parameters) comprise at least one of (A) time of year, (&) 
latitude and (C) altitude to generate an estimate of the zenith tropospheric delay. - 

As particularly described, the non-meteorological parameters are employed for 
determining a mapping function relating tropospheric delay at a given elevation angle 
25 to the zenith delay gen^-ated by this first model. A moping function using a 
continued fraction function based on the three non-meteorological parameters such as 
used hy Niell may be used for this process of generating the mapping function. 

The second model may be based on NWP or real time meteorological data, which is 

Currently the cell size of NWP data is 



10 



IS 



20 



25 



of the order of 12 km, but this is expected to reduce in me future allowing finer 
detenntnaiion of the model and mapping function. Using the NWP or other 
meteorological data (coUectivelyNWPO) a ray tracing technique is used to determine 
the paw of the satellite signal through the troposphere to the user and hence estimate 
the delay from a direct path. Knowing the delay, a second set of mapping function 
parameters may be derived by fitting the results to the format of the mapping function 
used in the first model. 

Recognising mat the wo models have generated two values for the same set of 
parameters, it is only necessary to communicate to the user the differences in the form 
of corrections to the four parameters estimated (zenith delay and me three parameters 
used to generate the mapping function) in me first model. Depending on the non- 
homogendtyofmetropospheieitmayonlybenecessaryto generate a subset of the 
parameters although the most important is the zenith delay followed in order by A, B 
and C. 

To be of use the data must be commimicated to the user within the time for which the 
NWPO is valid. As noted above the NWP data is defined over a geographic grid. If 
the standard deviation for the change in zenith delay for a cen in the grid in a period t 
equals the inaccuracy in the second set, then for a workable system (where the 
measurement period is less man changes in txopospheric conditions leading to 
misleading results) the data reduction needs to be effective to reduce the data to an 
amount which can be transmitted over the capacity of the data link in a time f 
Eubstantially less than, t, and preferably no more than t/4. 

The derived correction parameters are distributed according to location over at least 
part of the earth's surface and may be communicated to at least one geostationary or 
non-geostationary satellite or terrestrially based system. The GNSS satellites may 
also be used for mis communications link. To enable the use of low capacity data 
links it is preferred to compress the distributed data prior to transmission. Such 
reductions may conveniently be performed using lossy image compression (each as 
JPEG) or sub-sampling techniques. 




Further details and advantages of the invention will be evident upon a reading of the 
appended claims* and upon a reading af the following article. 




Tropospheric Delay Modelling and Correction 
Dissemination using Numerical Weather 
Prediction Fields 



Matthew Pow, James BuIGncr and John Owen 
Deface Science and Technology Laboratory 



ABSTRACT 

With foe modernisation of GPS and the introduction of 
Galileo, troposphere enore are co become an 
increasingly signrfiewt proportion of the GNSS user's 
error budget* This paper invwdgsues the spatial and 
temporal properties of meteorological features on the 
tropospheric anas experienced by GNSS use*. A 
regional leal-thnc tropospheric correction service could 
be provided flom lefiactrvity profiles, derived fiom a 
numerical weather prediction ttpML An analyst is 
described Of the bandwidth requirements ffr real-time 

twposjuteric correction disscmastwn ud of the twite- 
. eff between bacdwidih and accuracy. 

Whilst fl» dry atmosphere variability, both spatially and 
temporally is small* water vapour variability can he 
J«E»- High-resolution mesoscale numerical ^rather 
fields enable analysis of die abrupt spatial and temporal 
changes of atmospheric moisture associated with weather 
fronts. A tftree-dhnensional r*y-u*einfc technique has 
been developed to compute tropospheric: delays from 
aaesoscale numerical weather prediction fields for 
denned satellite elevation aod azimuth angles. 
Quantified examples of abrupt spatial, including 
axumidnl. asymmetry and temporal tropospheric delay 
variability are give* and statistical aoaryas of their 
likelihood presented. 

Once die tropospheric delay can be established, the issue 
of the data compression of regional tropospheric 
informatta* to form a practical service needs to be 
addressed. Diffineut strategics for domination of 
re^onal tmposphcrifl correct ions are discussed and die 
txade-oiT between accuracy and ccncchna service 
bandwidth described. 



1. WTRODTJCTION 

Tropospheric delay corrections can be derived from 3 
priori atmospheric mod^ » atn»spberie m 
real-time meteorological data or, gives sutSekut timev 
estimated from GNSS code and carrier observations. 



This paper investigates the potential of hhancrical 
Weather Prediction, NWP, models to provide ^ 
tropospheric corrections. 1 

The 3*D atmospheric profile mforrnarion resulting fiom 
the pf OCesa of numerical weather prcdicrioii offers the 
pcneatial of providing a regional or global tropospheric 
correct on service to GNSS users. La effect an accurate 
atmospheric refiactivity field, and therefore rtopospberic 
delay, is simply a by-product pf the ouoicrical tveaflier 
pmo^ctkni process. An accurate estimate of tropospoeac 
delay can be obtained by ray-tmcmg through the 
refractfvity field The accuracy of a NWr^erivcd 
tropospheric correction is almost eorirtly dependent upon 
the ability to taerrrrine atmospheric water vapour 
content 



2. 3D RAY TRACING 
1A Propagation Theory 

The speed of propagation of a* elcctrornagnchc wave 
through 4 medium can be expressed in terms of the 
refractive index; ft, defined tab* the ratio of the speed of 
light through ftee -space to the speed fiirou^ the incdr^ 
(50,3-1). 

c t>i) 
v 

where? 

n is the refractive index 

c is the speed of light ia fixe space 

v is the pftrpagation velocity 

The GNSS troposrAcric time delay, ignoring rclaiivistic 
erects, is defined to be the propagation time of the 
GNSS signal fiom the satellite to the user ramus the fiee 
space p ropaga t io n ti™: 
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where 



x is the distance along the propagation path. 
The first integral is along the curved 
propagation path; the second integral is along a 
geometric straight path. 



The differential equ^on" dc«?n*OJf^ thi curved ray paih 
can be ejgjrcsacd, m cartc^n coordinates^ as (3); 



v» 



(2-3) 



where r = ffa) is the vector describing the ray path, 
jt is the length of fie curved ray path up tp r t 
/t is the refiactivity scalar field. 
Vn t a vector field, is the gradient of ru 

The differential equation can be expanded as 

ay «^ ^ jfay * _ 

A firat order ordinary ififibttnoal equation (ODE) with 
known initial values can be solved using numerical 
methods: for example Rwage-Kutta ox Adams-Moulton 
methods. Higher order differential equations can be 
solved numerically by rewriting xhem as an equivalent 
system of first aider equations. Using the substitution r } 
' — r and =» V (die first derivative), the ray path 

differential equation (2-4} can be expressed as an 
equivalent system of two first order di£fcfexttial equations 
2-5 and 2-6: 



4r 



(2-5) 



(2-6) 



The detcminsrina of the ray path therefore amounts to 
the solution of a system of two ODEs with initial values. 
Matlab's inaction odedj was used to solve the problem: 
a Rungc-Kutta method with adaptive step control 
co nsistan l with user defined tolerances. 

The numerical weather prediction fields to be used are 
expressed in a spherical coordinate frame, it is 
computationally convemeur, therefore, to generate the 
refiwtiviiy gradient in spnerieal coordinates (r,$ff), 
which can be converted . into local curveilinear 
coordinates (u,v,w) using, the following transformation. 



i— u + 
dr 



SL W (2-7) 



A further rotational transformation is then applied to the 
local curveiliiicr coordinate frame to give the gradient in 
a feed cartesian feme (st.y,* in Fig. 2-1) suitable for 
numerically solving equations 2-5 and 2*6* 




Kg- 2^1 RcnaCuVity Geometry 



The tropospherie delay can >e computed as: 



where a, b and c are as shown in- Fig. £-2. Point b 
corresponds to the point at which ray curvature 
and reactivity is negligible, In this work 
assumed to be above an attitude of 7Gkm, 




Fig- 2-2 Ray Path Geometry 



2-2 Eskinple Propagation Paths 

An example propagation path for a satellite at an 
elevation of 5° is shown in Fig. 2-3. The co-ordinates 
used are and y</ as shown in Fig. 2-2. It is noted that 




the majority of the curvature or the ray-path, occurs^ 100- 
zOOfanftoTOthpww- It is the rotation of tbe reactivity 
Hiadirat vector field, caused by canh cujvanne, which is 
dominant cause of ray cwvaturei tamer than the rate 
of change of lefiacdviry scalar field. The apparent 
elevation angle, the F a P a S ation P 3 * ^ 
location, is approximately **■ a to eteVatV ° n 

angtefwm satellite to u?er of 5°. 

It is show, in Fie- 2-3. that hi extreme ionosphere (TEC 
=10OT does n°* result io a significant change in ray angle 
ra SSy into the trep«pibew therefore, the assiimpdcm of 
&e addition of (raposphcie io ionospheric delays is valid 
as ihe coupling effect is minimal 




Fig. 2»3 Example Ray Deviation at 5* 

Example delays and additional path tang* one to ny 
curvature are given in Fig. 2*4. The increase a padi 
length due to ewvwnw is approximately IScm for a 
satellite at 5° elevation. 




2.3 Equivalent Angle Approximation 

The method of dtterttmting the tropo^lwrie delay , by 
numerically solving a system of ordinary differentia! 
eqmnans, » ctoopwotionaUy expensive. A &»ter 
method is to integrate the refractive index |W»e* 
straight geometric path; however this jntegtaoon wffl 
result in an owiesthmtion of tcapaspfaene delay, xn 
reality me ray bends so as to mnrnn ise the propagation 
time (Fcnnar S principle). 

There exists an equivalent elevation angle, lying btfween 
the true (geometric) and apparent elevation angle, that 
results in the straight path integral being equal to the 
integral along the acta*! curbed path fiom satdiiw to 
user. The elevation correction to he added to the uue 
elevation to obtain the equivalent angle, varies win ww 
altitude, but is largely Unchanged vrifh atmosphere (rig. 
2-SL The standard deviation between the equivalent 
straight lias integral and me ray-traced TOpospherfe 
delays for devatinna between 3 and 10 was found to be 
rO.04%. 




Fig. 2-4 Example T>elay and Excess Curved Path 



1MOMHK PM0IMfl 

Fig. 2-5 CeirtCUOn UJ Teue Elevation Angle. 

3. NUMERICAL WEATHER YIELD PROCESSING 
3.1 Introduction 

Nunxrical Weather Prediction models mrecast the 
evolution of ain»spheric physical processes by applying 
governing equations, including die conservation of mass, 
jnomeinnm and energy. Tbree-dhnenainnal fields of 
coXuous variables including humidity, pressure, 
tenpeuume sod velocity are numerically grocq«d and 
meteorological features, hKlodmg «**ber fionu, am 
secondary derived propemes. A variety or 
nwasuremeats can be input mm me numerical rondel 
including sutfece, radiosonde and satellite observations. 
The water cyds is modelled including the effects of 
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terrain moisrure, sea surface rennpetfiturc. cloud 
fonnatfoa andpxejzjptf&iiOnu 

.• 

There are - a variety of eo-ordinaic systems in use for 
numerical weather prediction; they can be either based 
—upon- -grid of fixed points-or - spectral cwEfidems- - 
Physical height is rarely used as the vertical co-ordinate 
ia NWP: pressure; signa-pressure and hybrid co* 
ordinatea arc more convenient- Sigma-prcssure vertical 
co-ordinates define a level trj be a surface of constant 
pressure divided by surface pressure and thus have a 
terrain following property. Hybrid co-ordinates are 
retrain following at low levels* bus lend continuously to 
pressure surfaces at upper levels* 

Numerical models can be global or of limited area. 
Limited area hi gh-res alutioii models are often termed 
mesoscale models as they reflect mesoscale 
meteoiulogseal fisatmes: weather patterns of Jess than. 

lOOkm in size. 



3^-Mcso«ale Resrrfuiinn Model 



The rnesoacale data used is this study was obtained from 
the British Atmospheric; Data Centre whh permission 
from the United Kmgdcru Meteorological Office, 
UKMO. The UKMO mesoscale Unified Model uses a 
hybrid vertical co*ordinare sysrem with a spatial 
resolution of 12x12 ken. The model covers the area 
shown in Fig. 4-2. Specific humidity is modelled at 35 
hybrid levels; temperature is modelled at 38. To ray- 
trace through the numerical weather field the hybrid 
levels are convened into pliysica) height by assuming 
hydrostatic equilibrium. The impact of deviations from 
hydrostatic equilibriunt caused by mountainous terrain 
on lunar laser range measurements has been estimated to 
be less than 1cm for elcvaliod angles as low as 20 D [81 
die impact dec to coavccrrvc storms is not known. In 
assuming hydrostatic equilibrium, local values ibr 
gravity w used to reflect gravitational variation whh 
tanm& and altitude. 

Once the NWP fields are fixed in 3-D space, the data can 
be re-sampled, using a coznbinstion qf linear and log- 
linear interpolation, at any physical height. The 
atmosphere: was extrapolated from the top of the NWJP 
field (approximately 35km) go an altitude of 7Qfan 
Wffnng a- log*ljncar relationship between rc fiac t JYrty 
and altitude The atmosphere can then be tepro sentc d as 
a rtfraoivicy field, from which the refiraeriYity gradient ia , 
nutnencaxry comptuca* 
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Selectivity 



The atmospheric reiracuVtly can be derived from 
atmospheric pressure^ wafer vapour partial pressure and 
temperatu re fields. Atmospheric reJrsctxvity can be 
divided inco dry (hydrostatic) and wet components, A 



simple two-term expression, given in [2], has a 
lefcaetivity Accuracy of 0,5%: 



jVe|0*(ir-l) (3-2) 



where; 



tfistherejraetivity. 
P is atmospheric pressure (millibar) 
Eis the water vapour pressure (millibar) 
T is temperature (Kelvin) 



A review of more sophisticated three-term 
including gas ciinnncssfcjQity eiTccU, possessing greater 
accuracy, is given in [5], 

3.4 Future Numerical Weather Prediction 

- -The-rapid-improvement-m -power- and — 

advances in numerical techniques has seen the ECMwp 
model's hozizonlal spatial resolution reduce from 200km 
to. 26km in the period 1979-2001* This process » 
expected to continue and it is predicted thai by 2015 
lirnitcd^area model resolution could be ikm and the 
global model less than J 0km [ 1 1]. 



4. MESOSCALE METEOROLOGICAL FEATURES 

4 J Relevant Meteorological Features 

Meteorological features that possess a large spatial 
and/or temporal variation in tzopospheric delay will 
impact (he accuracy of NWP-derived tropospberic 
corrections and the bandwidth required ibr 
dwscroiuatfori. The temporal and spatial variation in 
hydrostatic refiaerivity b generally small, whereas 
. meteorological features associated with rapid changes in 
atmospheric moisture significantly impact . the 
aecuracy/bandwidrh relationship. 

The spatial and tcrnporal scales of certain weather 
systems are described in Fig. 4-1 (fiotti [!}). 

Meteorological features smaller than the resolution of die 
numerical p r e dic tion model will not be accucately 
reflected in the NWF -derived ox^cisj&efie correction. 
The 22x12 km resolution fields are able to reflect 
chaiecter^bca of weather fronts, however smaller 
features such as local storms and tornadoes would not be 
well represented spatially or temporally. The extent of 
: the weather Matures in altitude is critical in assessing 
their relevance to tropospherie delay anomalies. 
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Fig, 4-1 Spatial Scale and Lift Span of 
MeteorotoSKal System [1] 



4.2 Weather Fronts 

A w»*Cl Sent TO*5 the interfile between air: 
defined as a large body of air whose physical prnperues 
fen»ly umfcrra horizontally hiuntads oJ 
Stotnenk The front eaa raa* d* occurrence of abrupt 
chance in atmospheric nwistm^ tcn^^ ^ 
therefore refrachvity. firats can be dmded fato thro, 
classifications: warm, cold and occluded. 

Each ft™** zone fa imft l uc however tocy do pMScs 
some general characterise, frontal zones teal to have 
B width of a few lOkms. A w?rm front represent* 
jitojing edge of an advancing warm ™«. ™ 
BcnSy has a gentle slope (inclination^ WT»?»ny 0.5 - 
PT< =old ftooB, the leading «Jg« rf "dvanmng cooler 
air masses, are typically Alined » V. Another 
characteristic ofme fiont is its remove speed. CD ^ i frc ° te 

move much faster than warm firams, tip » SOantor 
compared to lOkmfer «^^J5^J^* ^ 
ftont catches np with a www ftw* * undcrctas tue 
warmer less dense air and the resultant, more complex, 
feature is tensed an occluded front Occluded fonts 
contain cbaraoerisdcs lhal belong to both coW and warm 
fronis. 

The most rapid change in tiopo^herie ^.j^ff » 
occur when sauttiu elevation and fiont mdinanon are 

equal. Generally. In ^W*^*?K2l22 
to 50 kilometre* per boor and can result m zenith delay 
! of 3cjWbonr (12J. 




Fig. 4-2 Synoptic Chart Weather From Feature 




Fig, 4-3 Meteosat Image 19VD1/ZDQ2 



An ocample of a weather Cental system shown m the 
^oST chart in fig. « (ootamcd torn to 
Krzentrak website). The ™^ U S? 
bighiianted Unified Model co^gc Hre^ mclud« 
3uded, add, warm and elevated coM fiontfc iw 
S^nTis associated with abrupt changes » 
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Fig. 5-1 Azimuthal Asjymmetry Feared Event 



atnwspberfc moisture content, which can be seen as as 
abrupt change in cloud coverage, Fig. 4-3. ft Mucosal 
image obtained, -from Dundee University; and specific 
humidity, Fig* 54* aNWPxsiesoscale image. 



5- AZIMCTHAL ASYMMETRY 

The accuracy of the- widespread assumption that she 
Atmosphere is symmetric in azimuth can be tested by 
examining [he variation in troposphere delay for a given 
elevation angle, A data set. cornpiising UKMO 
mesoscaie NWP fields is two-hourly iniervals for the i" 
and 15* of cadi month over a period of one year, was 
wed to generate azimuthal asymmetry amplitude 
statistics. Four user tactions were used dear the 
mesoscaie region centre; each location -was tested in 
in cremen ts of 10 degrees of attrnufli *t an altitude of 
SQOm above usen^sea-leveL Fig. 5-2 shows the 
percentage variation in tropnsphcric delay from the mean 
delay for ricvatwn angles of 2. S, 10 and 20°. The 
^trrryfrrrri, deviation of azimutbal asymmetry expressed as 
a percentage for each elevation is 0.21, 0.12, 0.07 and 
0-03% respectively. The amount of azinmtha] 
asymmetry is small in station to the current accuracy of 
the NWP-dcrrived utrpasphtnc corrections: a standard 
deviation of 0.4%. For users requiring the highest degree 
of accuracy, a u u r nu fbai asymmetry information would he 
of value, however tor most users, the cost of 
disseminating ™rrmrhnl asymmetry information 
outweighs the small gain in accuracy. As the accuracy of 
NWP-derived corrections improves, mainly due to 



improved model resciutian, sirfmnfrfflrf asymmetry is 
likely to become a more significant limiting factor 




Aztmufhal Asymmetry Pafconcaoo 
Fig- 5*2 Azzmmhal Asymmetry Histogram 

The distribution of the measure of asymmetry has long 
tails. Extreme levels of asymmetry occur when a 
weather front is overhead. Fig. 5-1 consists of a specific 
humidity horizontal slice at an altitude of 2km, a specific 
humidity vertical sKec and polar plot indicating the 
percentage of azimuthal asymmetry. The polar plot 
displays the percentage change in delay from the mean 
delay lor that elevation. The concentric rings are for 
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elev^on angles of 1, 2, 3. 4. S. 7» 10 and 15". TO 
variation of asymmetry caused by ample weather fionn 
can be accurately modelled as a sinusoid. Sinusoid 
parameters could than used by users requiring the highest 
accuracy of tropOSphttic corrections, as suggest in [7]. 

6. ELEVATION MAPPING 

In modelling troposphcric delay it is convenient to relate 
the troposphcric delay ac a given elevation angle to the 
jeoi* delay (4^) by means of a mapping function 



(6-1) 



The nydrosaiie zenith delay can be accurately modelled 
riven a surface pressure measurement [6J, howevet lie 
w« zenith delay can not be accurately detenmned fiom 
surface hnmidity measurements, as they aie not 

lepieseaiarlve of Ihe above atmn sphere. The expression 
for fins Uopospherie delay at a given elevation angle can 
• be defined as: 



It is noted flat care must taken when applying the 
simplification of die superposition of hydrostatic wet 
atmospheric delays: the propagation path is dependent 

Mpga both P '"""^ w*»f ertmnoilClltS. 



The w»tioa of troposphcric delay vtith elevation angle 
(Fig. 2-4) can be efficiently modelled by a continued 
fraction expansion. Niell hydrostatic and wet marbling 
fhncdons [10] powew a high degree of accuracy wiawuj 
the need for prior meteorological uafonmtron. Niell 
inapping factions use a continued fraciion function 
based on three pannneters a, b and. c, which arc * 
foiKrion of time of year, latitude and altitude. N«» 
mapphsg functions are likely to result in significant 
additional euros compared to the aceumcy of NWP- 
derived zenith correction*. A more accurate direct 
mapping fiuteuon can be used based upon NWP ray- 
traced results [4J, 

For the masons of dtocmnatwn efficiency, it was 
decided to we * rhree-WSin continued faction 
tppnsxkDaciott to u» total (we* and dry) mapping 
fraction determined from my tracing NWP field* 
(equation W), Parameters tto. *n *Od c a amount to an a 
priori NicU mapping fraction for the tool (hydwsiatic 
and wet) Uopospbeifc delay, «d «re a function of time of 
year, latitude and attitude. The pAtaxnelcTS An, M r and 
At? are a cdritctiuntepai to be broadcast. The correction 
parameters are determined from a fitting process such 
that the sum of die squares of the residuals between 
equation (6-3) and the ray-traced (truth) is inminased. 



Such an opproTdxnatioxi round, to be in agreement 
with the ray-traced results (typically «OJ% lesidual 
error}. 



1 + 



m(c s a 1 6,*) = 



(6-3) 



sin(c) + 



7. SPATIAL AND TEMPORAL VARIABILITY 

The rate of infonnawn requiied to accurately 
t fiqenmate troposphcric correction InfuiUMiion. m * 
regional basis. is dependent open the spatial and 
temporal variation of tropospherie delays. 

7.1 Temporal Resolution 

The impact of temporal resolution on *K accuracy of 
NWP -derived troposphcric corrections can be analysed 
hy comparing the change in zenith delay at time to and 
t+AL Using the data set as used for the aximuthal 
asymmetry analysis, ihe changes in zenith delay were 
computed for At values of L 2, 4 and 8 hours for aU 
points in the NWP coverage area. The standard 
deviation of temporal variation, expressed as a 
percentage £bt each &t value. * 0-475, 0-594. 0.79S and 
l.06S% respectively. The correspandiiig histo&rams are 
given in Fig, 7« I. 




pig_ 7_x Temporal Resolulian Histogram 

The impact of tetstnotal resolution is djgp^jj 
crfEripn to the current accuracy of the NWf^ 
rxqpospticrie cmrcctions: a standard deviation of 0-4%. 
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Fig- 7-2 Front Progression. 



A significant error additional error is therefore jroduced 
if ten^ora! resolution exceedsj hgmr. 

The distribution of the temporal variation has long talis. 
Fig. 7-2 displays the percentage change in zenith delay, 
far At values of 1, 2, and 4 hours, as a winter from 
traverses the UK. lac weather front has resulted in a 
tnaxSiamn rate of change of zenith delay of 4crn/hour 
(J .6%) measpred between V and to + 1 hour T 



7_Z Spatial Resolution 

The impact of spatial resolution on the accuracy of 
NWP-derived tropospberic eotrecdons can be analysed 
by comparing the change fn zenith delay for locations 
separated by fixed distances. . Using the data set as used 
. for the a zffittlthal ai^nimetxy analysis, the changes in 
zenilh del^y were computed for spatial separations of 1, 
2. 4 end S tnesoscale . grid points in North and East 
duecnons (equivalent .to 12, 24 f 48 and 9 go» 
separations). The standard deviations of spatial 
variation, expressed as a percentage for each spatial 
separation value, are Q.13B, 0.213, 0.352 and 0376% 
respectively.. The cotresponrfnis bgtogrgms an? given in 
Fig. 7-2, Spatial sesoJuttgna. larger than the icsoludon of 
die mesoscale model. 12 km, axe therefore found to result 
m significant addition TOpospheric error*. The spatial 
variation in troposphcric zeuith delay has a close 
relationship to temporal vacation. As a weather front 
passes over, for example at 40 km/hour in the east 
direction, the difference in zenith delay observed at e» 
and ta+l hour is similar to thai observed between two 
points separated by 40 km in the east/Vest direction. 




-1 -ftS -OS 42 ft OJ? 04 ojj, u I 

Zenith Delay Change Percentage * 
Fig. 7-3 Spatial Resolution Histogram 



& BAJSDWffJTH HEQUDUeirtENTS 
8.1 Quantisation 

To achieve ctemt rroposphcric correction 
dissemination, it was decided not to model €oe 
ataosphflrte, for example wfth pressure, humidity and 
lapse rate nafiminatEofl, but to supply the minimum 
information needed for accurate delay derctmimion. 
Similarly* for reasons of bandwidth efficiency* zenjih 
delay and mapping function parameters can be expressed 
as corrections to a priori models. The zenith delay, 
disseminated as a percentage cnrrrrtfiw to the RTCA 
model [9]. can be coded using & bits spanning the range 
-10% to +10% (consistent with [13]} with a resolution of 
0.0781% and therefore nrininaal ajxtatistiofl noise. 
Similarly, mapping function continued . fraction 
information can bo r ^ ; 1^rmiir' t f^ as a correction to a 
Niell-like total delay mapping function. It is proposed 
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una mapping function correction parameters can be 
adequacy described by M Bta (7. 5 and 4 for 
parameters £*r, Aft and te, Equation 6-3). 

&2 Image Comprefrinn 

To investigate the Impact of qwntiwtkm and 
compression an the accuracy of disawajiawa 
HoiwSBuWiC delay corrections, an example zenith delay 
J^Lved ftom UKMO tncsoscateTTWP fields, ™tii 
an active front was quantised and Con fess ed. T» 
accuracy of Che recovered (disseminated) correction was 
then compared to the original values: derived from ray- 
traced aumttieal weather prediciion fields. The resultant 
grayscale bitmap iruage, 146x142x8 Bits, and three 
JPEG images, with successively ukzc^c^k^ 
ratios, are shown in the first tour images of Fig. 8- 1. me 
size of the Wut bitmap image, * zipP*d (loss-less 
compression) bitmap, and flaw JPEQ images are shown 
in Table 6*1. 



25 



Type 


Size 
Bytes 


Comp. 
Ratio 


Noise 

fstd) 


Worst 
Value 


BMP 
Zipped WW 


28014 
11212 


1 

is 


0.023% 
0.023% 


0.05% 
0.05% 


JPEG90 
'"" JPEG30 


231S 


12.1 
24.6 


0.042% 
0.093% 


0.3% 
0-5% 


JPEG1Q 




38.2 


0.28% 


" 2,5% 






— 4 


0.035% 


035% 


H^t rat i ' ' 




16 


0.063% 


0.5% 



TflWee-l 



The bitmap <BMF) and m * 
quannsatmnnOiseofO.23%. ^/^^.^ 
compressed the bitmap image by a fccwr of 24.6. but the 
additional quantisation and compression noise is smau. 
less than 0.1%. compared » *e current accuracy of 
NWP-deoved correction (04%). It is noted however flat 
simple onder-sampling every fifth pixel, then using 
hifoear interpolation between, missing pixels, raaBsm a 
compression actor of 25 and less noise Am JPEGSft for 
this image JPEG image compression oflft* ™> advantage 
over a reduced sampling ate. Further compression 
resulted in substantial reduction in performance as can 
been seen in Fig. 8-1. 



Type 


Sw 
Byte* 


Ratio 


Note | 


Wont j 
Value 1 


BMP 




I 






Zipped BMP 


80404 


1,8 


0.025% i 




1 JPEG90 


ZQ1$5. J 


7 




JPEG50 


8457 


16.7 j 


0.17% 








4 


ai?% 


3.3% 




noise standard 
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Fig. 6X2 Example Zenith Delay Relative to HTCa Model (%) using UKMO Global Model 



An wwxnpJo UKMQ global NOT greyscal* bitmap 
image and two JPEG images, with successively 
increasing compression nrtws. 9bown in the last three 
images of Fig. 8-1- U» original iwnkb delay, expressed 
as a percentage change from the RTCA model, is shown 
in Fig. 8-2, The global model has a resolution of 5/6 and 
5/9 degrees in longitude send latitude respectively: 
producing a global image of 432x325x8 Bits. The 
JPEG90 image coicprcssw the bitmap image "by & factor 
or 7 (Tabic 8-2) with compression now of Jess than 
0,1%; tbc JPEG5Q image result? in significant additional 
zenith delay errors as can be seta in Fig. For the 
global image the suh-sampiwj technique, with a 
compression fcetor of 4, results in substantial additional 
compression errors. 

Assuming that mapping function parameter corrections 

can be similarly tomptess 

^* Bandwidth e.g. 209 BeVs global data disseminated 
in lOi 



on the error introduced due to 



9. SUMMARY CONCLUSIONS 

An accurate atmospheric mrmctivity profile,, ax*d 
therefore trnpospherio delay, is simply a by-product of 
the nuBjcncal weather jAevnctian process. A vector 
calculus based approach reduces the problem of 
propagation paifo detennjnatjora, through the numerical 
weather prediction derived refiactiyiiy field, to the 
solution of a system of ordinary differential equations 
suitable for solving by numerical methods. 

The accuracy of iropospheric corrections used by a 
satellite navigation user is dependent upon the accuracy . 
of the numerical -weather prediction (currently 0-4% 
standard deviation for the UKMO mcsoscale model) and 



of correction 



For the data set used in this study, errors due to the 
omission of azimuths) asymmetry effects were found to 
be smell in comparison with numerical weather 
prediction errors. Only for users requiring Ac highest 
degree of accuracy would wmnthal asymmetry 
information be of value. The impact of temporal 
resolution is significant in relation to the current 
accuracy of the NWF-derivcd tzoposphcric corrections. 
A significant additional error is therefore produced if 
temporal resolution exceeds 1 hour. Spatial resolutions, 
larger than The resolution of mc mesoscale model, 12 km, 
are therefore found to result In significant additional 

troposphetfc errors. 

For lessons of bandwidth efficiency, zenith delay and 
mapping function parameters cad be expressed as 
corrections to a priori models. 8 Bits con used to 
describe iropospheric zenith delay with minimal 
quantisation noise; a further 16 Bits msy be required to 
disseminate the elevation mapping function expre ssed as 
continued fraction correction parameters. 

The bopospherjc correction maps can be compressed 
with minimal reduction in image performance,, but a 
substantial decrease ha data size. A combination of 
efficient quantisation and concession techniques results 
in the data-rate for troposphere correction dissemination 
being suitable for transmission via Global Natfagtion 
Satellite Systems. 

The zenith delay and mapping fraction parameters can 
be used in assisting the process 
determination. 



14 



•- ** Bits ?ef Second results 



10L REFERENCES 

t Bany R. G- OwrfcylC JL~Al»M$totJt Wtatf>er*nd CTrnq*:"*, 
i mJ* v Hk Radto Refiaww Index: 1» Fcnw* 

3. a*™, mJW, e, taopw, sixth Mtica im 

4. ^^C^IS^^W. 3^ J. M» A. Hum, a. 

SS^Sw™ MvNnnfcer* Ft** 1593-1 «7. 

& Sort****** *t -Am»Vte«e OOroaio" 6* ^JW" 
nd atmtonhen a ndki of smiles, Tte Us* oT 

Artificial Mritffc* fa- Geodesy, Geopfays. Mono^ Sk* i* 2«7- 

T. a! Hcedo* * JU *««|teic «Smv^ 

asynsncCir on (fx* af jpoce xeedctis tUta", Jauraq4 W 

Gttttftptert Re»*Jl. VW. 102. B«, f>pa M^OT^C, 

* ttarcr. J, -a&ct* of Dwtofew Rem Hj*w«* 
BjuiEbrwrt on XtmbtpMt Corttttfertt W Satellite and I^mr 
Law Ranar M«wweai«n*\ Animal of Gsopftjpieal Roeareft, 
v«t 94. Hol B6. Fip» lUMWIfcV*** H*- 

* RTCa D0229C^ -MinntBim OperBiraal Ptafiroaac* Stt nflarti 
fbr dotal fosirtcmng Sysma ' Wid« Ana Aujpnennriau Sysvsa 
AiitencEquJjwtT^ NovaifcerXDOL. 

ia SMAE, -Gttbal Mg&it HatHa* ftp a* aumfram dffg 

■l mdio mbv^lii". Jnunal eT Orapfaywcal Research, Vol. 
tOl Jto. BZ, ftiB 3227-3244, fttewy 1996. 
u fichlr C -Ahiflfl terogical Weather ftqWw 300frgttft A 

UL c>WM5*T«Bte«ittQ^*n>c Eflfcct of WcaAcrJtoiatBiGre 
MeawremeflT, GTS Woett. tonovatf «v Ma* 1 993. 

13. Coffin*, ^ F- nd Unjky. K- B* ***** KoAial tropaspfcene 
prepajpfian delay: How lad cw IX PW6ttdm» «T ION 

O British OWW OWSlitt 3003U FtoMiso*! wi* peuiwooo , of 
Pftfto ce SCi<ttct awd Ttetnoiogy Labontoiy on behalT of CoaBOw 
HMSa 



CLAIMS 



-lv-- --A-methoto 

comprising the steps of generating a "first set of less accurate tropospheric delay values 
from a first model which is known per se, generating a second set of more accurate 
5 tropospheric delay values from a second model based on meteorological information, 
and developing a set of modified said parameters for use with said first model so that 
it then provides a set of delay values substantially approximating the second set, and 
expressing the modification to said first parameters as a set of correction parameters 
for communication to a user. 

—10 — 2: -A-mefhod-aocor^g^to^clmrfTl ^K&t^rm^^tmMel'is'^as^ im n<m-~ 

meteorological parameters. 

3. A method according to claim 2 wherein said non-meteorological parameters 
comprise one, or any two, or all three of time of year, latitude and altitude. 

4. A method according to any preceding claim wherein said first set comprises 
15 tropospheric zenith delays, 

. 5. A method according to claim 4 wherein the first model contains a mapping 
function relating tropospheric delay at a given elevation angle to the zenith based on 
methodology used by Neil to derive a parametric model of the tropospheric delay to 
satellites observed by the user. 

2G 6. A method according to any preceding claim wherein said second model is 
based on NWP data. 

7. A method according to any preceding claim wherein said second set is derived 
by a ray tracing technique. 

8. A method according to claim 7 wherein parameters generated in the first 
2 5 model are generated in the second model from the ray tracing data and the differences 

are determined for a geographical region. 
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9. A method according to ckiza 7 or claim 8 wherein a distribution of said 
correction parameters over at least part or all of the earth's surface is derived for 
communication to at least one geostationary or non-geostationary satellite. 

10. A method according to claim 9 wherein said distribution is subject to data 
5 reduction to enable a low capacity communications channel or data link to be used. 

11. A method according to claim 10 where said data reduction is performed using 
lossy image compression or sub-sampling techniques. 

12. A method according to claim 10 or claim 11 and including the step of 
transmitting the reduced data over a communications channel or data link. 

10 13. A method according to claim 12 wherein the amount of data reduction is 
sufficient to permit said transmission within a time substantially lower than the 
validity time of the said meteorological information, 

14. A method according to claim 12 Wherein the NWP data defines a cell size, and 
if the standard deviation for (he change in zenith delay for a cell over a period t equals 

15 die inaccuracy in the said second set, then the data reduction is effective to reduce the 
data to an amount which can be transmitted over a bandwidth of 200 bit per second in , 
a time tl substantially less than , 

15. A method according to claim 12 wherein V is no more than t/4. 

16. A method of any one of claims 12 to 15 wherein said communications 
2 0 channel or data link is to a geostationary or non-geostationary satellite. 

15. A method according to claim 14 wherein said geostationary or non- 
geostationary satellite retransmits at least part of said distribution to a user. 

1 6. Apparatus of obt aining data for use in a satellite positioning system or GNSS 
comprising first generating means for generating a first set of less accurate 
25 tropospheric delay values from a first model which is known per se and which is 
based on non-meteorological parameters,, second generating means for generating a 
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second set of more accurate tropospheric delay values fiom a second model based on 
meteorological information, and developing means for developing a set of modified 

" safd parameters fbFurewim saTd"fira moa^so thai nTto"^ ^ Via of delay 
values substantially approximating the second set, and said developing means being 
arranged to express the modification to said first parameters as a set of correction 

parameters. 

17. Apparatus according to claim 16 and including means for compressing said set 
of correction factors. 

18. Apparatus according to claim 16 or claim 17 and comprising transmission 
_means forjransinitting. said.setor- said-compressed- set-of-correctionfactorsto-a-uscr - 
via a geostationary or non-geostationary satellite. 
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